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ABSTRACT 

We examine the constraints on soft X-ray photon emissions from the reionization era. It is generally 
assumed that the Universe was reionized by ultraviolet photons radiated from massive stars. However, 
it has been argued that X-ray photons associated with the death of these stars would have contributed 
~ 10% to the total number of ionizations via several channels. The parameter space for a significant 
component of cosmological reionization to be sourced by X-rays is limited by a few observations. 
We revisit the unresolved soft X-ray background constraint on high-redshift X-ray production and 
show that soft X-ray background measurements significantly limit the contribution to reionization 
from several potential sources: X-rays from X-ray binaries, from Compton scattering off supernovae- 
accelerated electrons, and from the annihilation of dark matter particles. We discuss the additional 
limits on high-redshift X-ray photon production from (1) z ~ 3 measurements of metal absorption 
lines in quasar spectra, (2) the consensus that helium reionization was ending at z « 3, and (3) 
measurements of the intergalactic medium's thermal history. We show that observations of z ~ 3 
metal lines allow little room for extra coeval soft X-ray emission from a nonstandard X-ray sources. 
In addition, we show that the late reionization of helium makes it quite difficult to also ionize the 
hydrogen at z > 6 with a single source population (such as quasars) and that it likely requires the 
spectrum of ionizing emissions to soften with increasing redshift. However, we find that it is difficult 
to constrain an X-ray contribution to reionization from the intergalactic temperature history. We 
show that the intergalactic gas would have been heated to a narrower range of temperatures than 
is typically assumed at reionization, 2 — 3 x 10 4 K, with this temperature depending weakly on the 
ionizing sources' spectra. 

Subject headings: cosmology: theory — large-scale structure of universe — quasars: absorption lines 
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1. INTRODUCTION 

Extreme ultraviolet (EUV) radiation from stars 
within diminutive galaxies is the leading candidate 
for w hat sourced hydrogen reionization at z ~ 10 
(e.g., IWvithe fc Loebl l200l . This scenario has been 
the focus of the majority of t heoretical and numeri 
cal work on reionization (e.g., iFurlanetto et al 
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Finlator et al.ll201lD . However, other models may still 
be viable, such as models with a significant contribu- 
tion to reionization from active galactic nuclei (AGN; 
iVolonteri fc Gnedinl [2001 ) . The bright end of the AGN 
luminosity fun ction declines dr astically between z = 3 
and z = 6 (!Fan et al.l 120011) . precluding rare AGN 
(i.e., quasars) as the primary driver of reionization. 
Nevertheless , reionization by fainter AGN may still 
be possible (ISiana et al.l 120081 : iShankar fc Mathurl 120071 
iGlikman et al.ll201lL although see lWillott et alJ^OloTf ). 
In addition, at least a fraction of reionization could 
have been sourced by soft X-rays produced from inverse 
Compton scattering of cosmic microwave background 
(CMB ) photons off supernova-accelerated electrons (jOhl 
120011 ). by em ission produced directly in supernova rem- 
nant shocks dJohnson fc Khochfarl|2011[). bv high-mas s 
X-ray binaries (jPower et al.l 120091 IMirabel et al.l [20~lTh . 
by mini-qua sars dMadau et all I2004D ■ by ha lo accre- 
tion shocks (jDopita et a ll 1201 It IWvithe et all 120111). o r 
by dark matter annihilations (jBelikov fc Hooper! 12009). 

1 Department of Astronomy, University of California, Berkeley, 
CA 94720, USA; mmcquinn@berkeley.edu 



Even if stars reioni ze d the bulk of the inte r galac tic hy- 
drogen. lObl jt200lh. Uohnson fc Khochfarl (j20lTt ). and 
IMirabel et all (|2011l) argued that X-rays associated with 
the remnants of these stars should contribute ~ 10% of 
the ionizations via any of a few mechanisms. 

X-ray reionization scenarios have recently received re- 
newed attention (e.g.,[Haiman 2011), particularly for two 
reasons. Firstly, soft X-ray photons would have had lit- 
tle difficulty in escaping from z ~ 10 galaxies and ioniz- 
ing the intergalactic medium (IGM). This is in contrast 
to EUV photons, for which it i s unclear whe ther they 
would have escaped at all (e.g., iGnedinl 120081 ) . Empiri- 
cally, observations of z ~ 3 — 4 Lycc emitters suggest that 
a non- negligib le fraction of ionizi ng photons was escaping 
into the IGM (jNestor et alll201lD . Curiously, for Lyman- 
break ga laxies the escape frac tion is constrained to be 
smaller. iVanzella et all (|2010f ) placed an upper bound 
on the escape fraction of 3.4 < z < 4.5 lyman breaks of 
< 5 — 20%. Secondly, there is evidence that high-mass X- 
ray binaries (HMXBs) and ultra-lumino us X-ray sources 
(ULXs) anti-correlate w ith metallicity (|Crowther et all 
120101 iKaaret et aT]|201lD . This anti-correlation suggests 
that these X-ray sources - which already dominate the 
X-ray production of low-redshift star-for ming galaxies — 
were even more prevalent at high redshifts ( Mi rabel et all 
l20ll . 

A full or partial reionization by soft X-rays would have 
resulted in a different morphology of intergalactic H 11 
regions during this process compared to a full reioniza- 
tion by EUV photons alone. The mean free path of soft 
X-ray photons is much longer than EUV photons owing 
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to the strong energy dependence of the photoionization 
cross section. A photon with energy E 1 has a mean free 
path, Ahi, to be absorbed by hydrogen of 
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where xm is the fraction of hydrogen that is neutral. 2 
If a significant fraction of the IGM were reionized by 
E 1 > 200 eV photons, the inhomogeneous structure on 
scales of ~ 10 comoving Mpc that is anticipated by most 
models of this process would have been largely erased. A 
more homogeneous reionization scenario than that pre- 
ferred by models may be more consistent with limits on 
the k inetic Sunyaev-Zeldovich anisotropy from r eioniza- 
tion (jReichardt et al.ll2012t iMesinger et al.ll2012f ). High- 
redshift X-ray production is also important for the re- 
heating of the Universe and determines whether red- 
shifted 21cm radiation appears in emission or absor ption 
(jFurlanetto et al.|[200l iMcQuinn fc O'Leary|[20T2l ). 

This study investigates the constraints from current 
data on high-redshift X-ray emissions. Because the Uni- 
verse is optically thin to > 1 keV photons from z ~ 10, 
such photons produced during the reionization era would 
have streamed freely and altered the ionization state of 
intergalactic metals at z ~ 3 as observed in quasar ab- 
sorption spectra. They also wou ld contribute to the 
observed soft X-ray b ackground (jDiikstra et al.1 120041 : 
iSalvaterra et al.1 120071 ). In addition, the hardness of 
the high-redshift EUV and X-ray emissions determines 
whether the electrons of helium were ionized concur- 
rently with those of the hydrogen or well after. An 
early reionization of the second electron of helium ap- 
pears to be in conflict with the mounting evidence that 
helium reionization ended at z rs 2.7 (I Theuns et ahi 
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2002 ; iHui fc Haimanl 120031: iFurlanetto fc D ixon 
McQuinnl 120091; iShull et al l 120101; iBecker et al 
Worseck et al.ll201ll) . Finally, measurements of the ther- 
mal state of the IGM h ave been extended t o 4 < z < 6.5 
in th e past few years (Bec ker et al.l 1201 It iBolton et al.1 
120111) . redshifts more sensitive to the heating during hy- 
drogen reionization. Reionization by a harder spectrum 
than stars would have injected more heat into the IGM, 
potentially resulting in higher temperatures. 

Section [2] summarizes the leading candidates for sourc- 
ing the production of high-redshift soft X-ray and EUV 
photons. Section [3] discusses the constraints on high- 
redshift X-ray production from the unresolved soft X-ray 
background, from z ~ 2.5 metal absorption line observa- 
tions, from our knowledge of the reionization history of 
the hydrogen and helium, and, finally, from constraints 
on the thermal history of the IGM. In our calculations, 
we assume a flat ACDM cosmological model when nec- 
essary with n m = 0.27, h = 0.71, ct 8 = 0.8, n s = 0.96, 
and £lh = 0.046, c onsistent with recent measurements 
(jLarson et al.ll2011l) . 

2. POTENTIAL X-RAY SOURCES AT HIGH- REDSHIFTS 

2 The mean free path to be absorbed by He II and He I is 
comparable, with Ariel = 4 a:JJ eI cMpc and AHell = 5 ^JJell cMpc, 
assuming _E 7 = 200 eV and z = 9. Ariel has a softer scaling with 
&y than Ahi such that the He I would be the first species to be 
ionized by a background of just soft X-ray photons. 



Here we provide a brief synopsis of source populations 
that studies have argued could dominate the production 
of X-rays at high-redshifts and ionize > 10% of the Uni- 
verse. 

AGN: The abundance of the brightest AGN i s observed 
to fal l off rapidly with redshift at z > 3 (|Fan et al.l 
I2001D . However, the abundance of AGN that lie on 
the faint end of the luminosity function at z > 3 is 
not as well constrained and could conceivably have 
grown with increasing redshift (although see §2.1[) . 
The average spectral index of the intensity per unit 
frequency for low-redshift Type I qua sars is mea- 
sured to be a = —1.7 in the X-ray (|Tozzi et al.l 
2006) and is found to have a similar index in the 
far UV dTelfer et al.l(2002T) . 3 

The end products of stellar evolution: Previous 

studies have suggested several possibilities that fit 
within this class: 

1. from HMXBs - Emission from accretion onto 
the primary - a black hole or neutron star 
- dominates the X-ray output of low-redshift 
star forming galaxies. iMirabel et"aT1 (|2011f l 
argued that these sources are likely to be more 
abundant per unit of star formation at high 
redshifts and, thus, to have contributed to the 
reionization of the Universe. The population- 
averaged spectral index for HMXBs has been 
measured to be a ~ —0.7 at 1 — 10 keV 
(jSwartz et al.ll2004t iRephaeli et al.lll995f ). but 
it is unknown whether this power-law contin- 
ues to lower energies. 

2. from the cooling of supernova-accelerated elec- 
trons - At high redshifts, electronic cos- 
mic rays that were accelerated in supernovae 
would have cooled primarily b y inverse C om - 
ton sc attering off the CMB |On1[200l . [O 
(2001) argued that this mechanism could nat- 
urally yield 10% of the ionizing photons pro- 
duced by a stellar population if / esc w 0.01, 
where / OS c is the fraction of 1 Ry photons that 
can escape from within galaxies. This mecha- 
nism generates a spectrum with spectral index 
in Ie of —7/2, where 7 is the input power-law 
of the electron spectrum. It is measured to be 
w 2 for the injected cosmic ray spectru m in 
the Milky Way halo (|Strong et al.ll2004D 
subsequently assume 7 = 2. 4 



We 



3 If just 1% of the energy to grow into a black hole of mass 
10 5 Mq goes into ionizations and there is 1 black hole for every 10 
cubic cMpc of volume, black holes would have been able to reionize 
the Universe. (High-redshift quasar luminosity function measure- 
ments are only sensitive to AGN with > 10 7 Mq ; IWillott et al.l 
2010a). This mass density in massive black holes is a factor of 25 
times smaller than the SMBH density at 2 = IIYu fc Tremainel 
2002) and is consi s tent w ith the conservative theoretical estimates 
in lVolo ntcri ct al. (2008). In addition, if there were a relation that 
tied black hole mass to the total stellar mass in a galaxy at high 
redshifts (similar to the Magorrian relationship at lower redshifts), 
ILidz et al.l ((2007) showed that black holes would have produced a 
comparable time-integrated number of ionizing photons to stars. 

4 The fiducial values for X-ray production in Oh (2001) arc likely 
optimistic. Estimates for the fraction of the ejecta energy that goes 
into accelerated electrons find 0.1 — 1% rather than the 10% value 
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3. from shocks in s u perno va remnants 
Joh nson fc Khochfarl (|201lD estimated that 
direct emission from such shocks produces 
~ 10% of the ionizing photons that originate 
from a stellar population. Using models of 
shock emission, we find that this process 
results in a time-integrated spectrum that 
can be approximated with a spectral index of 
a ~ — 1.7 at < 1 keV and with a falloff above 
this energy. 5 

Dark matter annihilations: The electronic by- 
products cool by inverse Compton (or synchrotron) 
emission. This cooling results in a steady-state 
e+e- spectrum with equal energy distributed per 
log energy up to energies of order the particle 
mass and an inverse Com pton spectrum in Ir with 
spect ral index —1/2 jM cQuinn fe Zaldarriagal 
l20ll . iBelikov fe Hooped (|2009f ) found that "a 
generic weakly interacting massive particle could 
result in an ionized fraction of 10~ 2 and that 
recently popular models with an enhanced anni- 
hilation cross section could even result in a full 
rcionization. 

Some of the aforementioned sources could be partially 
obscured by a column of hydrogen and helium gas, which 
would further harden their ionizing emission. 6 

2.1. Direct Observations of UV and X-ray Sources 

Measurements of the quasar optical luminosity func- 
tion as well as of X-ray source counts place constraints 
on the abundance and evolution of the aforementioned 
sources. The standard lore is that the quasar lumi- 
nosity function falls off quickly towards high redshifts 
such that quasars were a subdomin ant contribution to 
the reionization of the IGM (e.g., iMadau et all 119991 : 
iFaucher-Giguere et al.ll200"8al ). However, this conclusion 
has stemmed from observations of bright, L > L+ quasars 
at z > 3. 

Recent studies have begun to constrain t he faint-end 
slope at z > 3 (IDiikstra fc Wyithd 120061: iSiana et all 
2001 IShankar fc Mathu rl 120071: IWillott et all 127)101 



Glik man et al. 1 1201 Hi. Of note, IWillott et al 



placed an order of magnitude more stringent constraint 
on the abundance of faint z ~ 6 quasars than other 
analyses, using the Canada-France High-z Quasar Sur- 
vey. With an absolute magnitude limit of AfusnA = —22 
plus ex trapolations to fainter magnitudes, IWillott et all 
(|2010H ) found that quasars could contribute no more 
than 10% to the ionizations required to maintain reion- 
ization at z = 6. If this result is confirmed and there 
is not a significant turn-up in the quasar abundance to- 
wards fainter luminosities, quasars alone cannot reionize 



used in this study (e.g., Kobayashi et al. 2004; Thompson et al. 

i2ooa v 

5 This calculation uses the spectral libraries from supernovae 
shock calculations of Allen et al. (2008), averaging over shock ve- 
locities of < 1000 km s —1 as outlined in Johnson & Khochfar 

poTm . 

6 A hydrogen column of TVhi has a photoionization optical depth 
of unity to photons with E < 1 (N m / 10 23 cm~ 2 ) 035 ke V. Signifi- 
cant obscu ration impacts 75% o f AGN (Tozzi et al. 2006) and many 
HMXBs l ILutovinov et al.ll25o5l ). 



the Universe. 7 In this paper, we will remain agnostic as 
to whether AGN reionization scenarios are allowed and 
investigate whether they can also be constrained with 
other techniques. 

At the faintest fluxes measured with the Chandra 
Space Telescope, star-forming galaxies become an impor- 
tant contribution to the ob served X-ray source counts 
(|Hickox fc Mar kevitchl 120071 ). In fact, if the faint end 
of the quasar luminosity function does indeed decline 
with increasing redshift at z > 3, the dominant high- 
rcdshift X-ray sources in the Universe were likely star- 
forming galaxies. There are some constraints on how 
the X-ray em i ssions fro m galaxies evolve with r edshift. 
INandra et all (|200l and lReddv fc Steidell ((2001 found 
that stacked z = 1.5 — 3 Lyman-break galaxies are 
consistent with the hard X-ray lumin osity-star forma- 
tion ra te relationsh i p foun d at z = (jReddv fc Steidell 
120041 ) . iCowie et al.l (|2012l ) reached a similar conclusion 
for galaxies at redshifts of 1 < z < 6. They found a 
fa ctor of 3 higher nor malization than the z = relation 
of lMineo et all (|2012D at all redshifts (but s imilar to the 
z = determination in IDiikstra et alll2012l ). with a fac- 
tor of few uncertainty on their determination at z > 4. It 
is worth noting that these constraints are on star forma- 
tion in the most luminous high-redshift galaxies, which 
may not be representative of most coeval star formation. 

3. INDIRECT CONSTRAINTS ON X-RAY SOURCE MODELS 

The following investigates additional constraints on 
these X-ray source models from measurements sensitive 
to the total amount of emission (in contrast to luminos- 
ity function measurements, which often probe only the 
brightest contributors) . 

3.1. The Soft X-ray Background 

This section revisits the constraints on X-ray reion- 
ization scenarios from observations of the soft X-ray 
ba ckground (SXB ) Such constraints were first discussed 
in IDiikstra et al.l (|2004l ). and an aim of this section 
is to update these calculations. The soft X-ray back- 
ground is one of the most constraining backgrounds 
for z ~ 10 sources because it is not altered by in- 
tervening absorption, in contrast to optical and ultra- 
violet backgrounds. Our focus is on the constraints 
from reionization scenarios in which only a few ioniz- 
ing photons produced per hydrogen atom were required 
to complete reionization, a differe nt region of pa r amete r 
space than originally studied in IDiikstra et all {2004). 
Both observations of the ~ 1 Ry emissivity of sources 
from the Lya forest at lower redshifts and investiga- 
tions of the expected number of recombinations dur- 
ing rei onization suggest that reionization was in this 
regime (lMiralda-Escudel200l iBolton fc Haehneltl l2007bl: 
iPawlik et al.l 120091 : iMcQuinn et al.ll2011l ). 

The SXB can be used to constrain the comoving num- 
ber density of ionizing photons, rii, emitted before the 
end of reionization. If this number density is greater 
than the number of hydrogen atoms, then there were 

7 We note that all estimates for the quasar ionizing photon emis- 
sivity use a single template to extrapolate from the source-frame 
far UV to the source-frame EUV. This extrapolation would be in 
error if spectral hardness anticorrelates wit h ultraviolet luminosity, 
as has been suggested in Scott et al. 20041 . 
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enough photons to reionize the Universe. In particular, 
in the absence of foreground emission, is related to the 
intensity of the z = background per unit energy, Ie, 
via 



47T 



c/(l+Z E ) 



13.6eV/(l+2 B ) 



dE-^exp[-T eS (E,z E )], (2) 

ill 



where exp[— r c g (E, ze)\ is the average transmission for a 
photon that travels from redshift ze, reaching the Earth 
with energy E, and here we have assumed for simplicity a 
single redshift z E for the emission. Our subsequent calcu- 
lations will assume E max = 1.8 [(1 + z)/15)] ' 5 a^ 3 keV, 
where xm is the H I fraction, which we take to be 0.5. 
This corresponds to the photon energy that has a pho- 
toionization optical depth equal to unity to travel a Hub- 
ble distance at this neutral fraction and redshift. How- 
ever, our results depend weakly on E mayi . 

X-ray photons originating at z ~ 10 with energy 
> 1 keV have a negligible probability to be absorbed 
(i.e., exp[— r c ff] = 1). Thus, since I E is measured in 
the X-ray, if we make assumptions for how Ie scales 
to lower E, we can constrain m. We adopt 3.5 x 
10~ 13 erg cm~ 2 s _1 deg -2 for the unresolved soft X-ray 
background int ensity between 0.5 an d 2 k eV. This value 
was d erived in iMoretti et al.l (|2003[ ) and iDiikstra et al.l 
( 2004) by combining both widefield and pencil-beam X- 
ray surveys and also by subtracting off the small con- 
tribution owing to Thomson scattered radiation from 
known sources. This ba ckground estimate is consis- 
tent with the estimate of IHickox k Markevitchl (|2007f ) 
from the Chandra deep fields, which was 3.5 (±1.4) x 
10~ 13 erg cm~ 2 s" 1 deg~ 2 between 1 and 2 keV. (We are 
assuming that all of the resolved X-ray sources are from 
after reionization. This is likely the case. In fact, the 
IHickox k~ Markevitchl 120071 measurement masked point 
sources detected with the Hubble Space Telescope in the 
Z-band, which incidentally would have missed masking 
sources at z > 6.5 owing to H I Lyman-series absorp- 
tion.) 

For a spectral index of a for the sources' intensity per 
unit frequency, the number of ionizing photons produced 
per baryon Ni(a) = rii/nn if all of the unresolved back- 
ground owes to reionization era sources is 

JVi({-0.5, -1, -1.5}) « {0.05, 0.8, 14} 7 3 .5 (3) 

{-0.5,-1,-1.5} 

/ ' T . /■; >, 

X 



1 + Ze 

10 



where 13.5 is the 0.5 



background in units of 3.5 x 10 13 erg cm 2 s 



2 keV intensity of the soft X-ray 

1 de s~ 2 ' 

and ze is the emission redshift. Even though we have as- 
sumed that all the photons were emitted at a single red- 
shift, the results are rather weakly sensitive to this choice 
for reionization era sources. Equation demonstrates 
that the SXB can rule out X-ray reionization models that 
satisfy a > — 1 as these have JVj < 1. 

Figure [1] illustrates in more detail the upper bound 
from the unresolved SXB on Ni as a function of a. The 
solid curves (dot-dashed curve) are for the case in which 
all of the photons were produced at z = 10 (6), assuming 
no absorptions from within the progenitor galaxy. The 
thin solid curve does not include secondary ionizations, 
while the three thick curves include them (which is the 
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Fig. 1. — Constraints from observations of the unresolved soft X- 
ray background (SXB) intensity on the ionizing photon production 
of high-rcdshift sources. Shown is the number of such photons 
per hydrogen atom a source population could have produced as a 
function of its effective spectral index in Ig between 1 Ry and the 
hard X-ray band, denoted as a in the text. The solid curves (dot- 
dashed curve) are the contours that saturates the unresolved SXB 
for the case where the photons were produced instantaneously at 
z = 10 (6). The thin solid curve omits secondary ionizations, while 
the the other curves account for them as described in the text. For 
the scenario represented by each curve, the parameter space above 
its curve is excluded. The "GRB DLAs" curve is the same as the 
thick solid curve but the photons are assumed to have been filtered 
by the host-galaxy H I column density distribution inferred from 
GRB afterglows (which models interstellar absorption for galactic 
sources). The labeled vertical lines show the fraction of ionizations 
that could result from the different source populations discussed 
in Ej2j using empirically motivated estimates for their power-law 
index. 



more realistic case). Secondary ionizations, which were 
ignored in the previous discussion, are included by as- 
suming that the ionizations occurred in a medium that 
was initially neutral, ignoring reco mbinations, and usin g 
the secondary ionization rates in iRicotti et al.l {2002). 
(Even during a patchy reionization scenario, X-ray ion- 
izations typically would have occurred outside of the H 11 
regions that surround EUV sources where the IGM was 
largely neutral.) This approximation allows us to map 
the previous estimates for m (eqn. [3|) to a boosted n% 
that includes secondaries. The differences between the 
thick and thin solid curves demonstrate that secondary 
ionizations only impact Ni for C 1. Otherwise, the 
excess energy per ionization is channeled primarily into 
heating the IGM rather than into secondary ionizations 
(jShull k van Steenberd[l985l ). 

The dashed curve labeled "GRB DLAs" in Figure Q] 
is the same as the solid curve but where the emitted 
ionizing photons are filtered by the host-galaxy H 1 col- 
umn density distribution inferred from gamma ray burst 
(GRB) afterglow spectra. It is plausible that this distri- 
bution is the same as the column-density distribution to 
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escap e galaxies from star forming regions (|Gnedin et al.1 
120081) . This GRB column-densi ty distribution res ults in 
/esc « 0.02 for 1 Ry photons (jChen et al.l [2007h . The 
labeled vertical lines in Figure [1] represent the spectral 
index expected for the sources discussed in $2] The 
constraints o n AGN are weak , whic h agrees with the 
conclusion of iSalvaterra et ail (|2007l ) that AGN reion- 
ization scenarios can be consistent with the unresolved 
SXB. We have conservatively assumed that the average 
spectral index of AGN is —1.7 as found for quasars at 
lower redshifts; either obscuration or a hardening of the 
spectrum with redshift (potentially stemming from the 
emission originating from lower mass central black holes) 
would tighten these bounds. The constraints from in- 
verse Compton from supernova-accelerated electrons are 
more significant. This mechanism could not fully reion- 
ize the H I and is limited to ionizing < 30% of the Uni- 
verse in the case with GRB damped Lyman- a systems 
(DLAs) . The same conclusion holds for any process that 
results in a steady state number of electrons per unity en- 
ergy with a power-law index of —3 in energy or harder. 
The steady-state electron distribution would be —2 for 
the electronic byproducts of dark matter annihilations in 
standard weakly interacting massive dark matter parti- 
cle models (Section [2]) . In this case, the constraints on 
Ni are even tighter (iV, < 10%); however, in this case 
the approximation of the sources turning on and off at 
a single redshift cannot hold, and the associated ioniz- 
ing emission is likely to also occur at lower redshifts. 
Metal absorption line measurements rule out a signifi- 
cant contribution from X-ray sources that are still active 
at z = 3 with a > — 1 ( t!3.2[) . If we use the estimate 
for the spectral index for the galaxy-av eraged emission 
of HMXBs and from UL Xs of a = -0.7 (Rc phaeli et al.l 
119951: iSwartz et al.ll200l . then HMXBs/ULXs can also 
only contribute < 10 — 20% of the ionizations. If their 
spectrum were harder than this index, as is likely at lower 
frequencies (especially since this emission is largely ther- 
mal radiation from a truncated accretion disk with char- 
acteristic temperatures of T ~ 1 — 10 keV), then the limit 
on their contribution would have been even tighter. 

We note that all of the aforementioned constraints im- 
prove proportionally to the measured unresolved SXB 
intensity. It is quite likely that a large fraction of the 
unresolved X-ray background was produced after reion- 
ization by galaxies below the current flux threshold of the 
observations, by the local bubble, or the warm hot inter- 
galactic medium. In future missions, lower frequencies 
(that are closer to the EUV band) would be interesting 
to target with SXB measurements, as we estimate that 
the foreground extragalactic H I absorption for photon 
energies > 200 eV has effective optical depth of less than 
unity (and the optical depth contributed by the Milky 
Way is also less than unity at > 200 eV). 

Figure [2] considers the type of evolution that is nec- 
essary for galactic X-ray emission (probably associated 
with HMXBs) to saturate the unresolved SXB inten- 
sity. 8 The curves represent the 2 — 8 keV luminosity 
per unit star formation rate, cx, needed to saturate the 

8 The calcu l ations shown in this figure are similar to those in 
Diikstra et al. (2012) except that we do not tie our constraints to 
the observed psfr(z), which suffers from incompleteness at high 
redshifts. 
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Fig. 2. — The 2 — 8 keV luminosity produced per unit star for- 
mation rate, cx, required to saturate the unresolved SXB, for dif- 
ferent cosmic star formation rate densities, Psfr- c x is given in 
units of erg s — 1 /(Mq yr — 1 cMpc - 3 ). The curves show the val- 
ues which saturate the unresolved SXB under the assumption that 
the sources radiated with a constant emissivity from z = oo to 
either 2 = 3 (solid curves) or to z = 8 (dashed curves). The thick 
curves are for a = — 1, and the thin are for a = 0. The dotted 
horizontal lines denote the star formation rate density required to 
maintain reionization at z = 6 for / csc = 0.1 (lower horizontal 
line) and / eS c = 0.01 (upper horizontal line; see text for details). 
The shaded region is the constraint on the relationship between 
X-ray luminosity a nd star formation rate from z = observations 
IMineo et al. 2012). Its width is set by the observed scatter about 
the mean relation. 

SXB for different cosmic star formation rate densities, 
Psfr- The curves assume that the sources radiate at 
a constant emissivity from z — oo — but the depen- 
dence on the maximum redshift is weak - to z = 3 
(solid curve) and z — 8 (dashed curve). The thick 
curves are for a = —1, and the thin for a = 0. The 
lower and upper horizontal lines denote the star forma- 
tion rate (SFR) for a Salpeter initial mass function re- 
quired to maintain reionization (balance recombinations) 
for a clumping factor of 5 at z — 7 for / csc = 0.1 
(requiring psfr ~ 5 x 10~ 2 M Q yr _1 cMpc~ 3 ) and 
f Rsr , = 0.01 (requirin g psfr ~ 5 x 10 _1 M Q yr _1 cMpc -3 ) 
(jMadau et all 119991 ). Estimates from the Hubble Ul- 
tra Deep Field fall on the low side of these psfr, fa- 
voring psfr ~ 10~ 2 M Q yr^ 1 cMpc -3 at z = 8 
and an order of ma gnitude larger value at z = 3 
l|Bouwens et al.l [20111 ). but such estimates miss fainter 
galaxies which likely dominate ps fr, especially at z > 6 
(jKuhlen fc Faucher-Gig ucrc 2012j ). The shaded region is 
the constraint from IMineo et al.l (|2012f ) on cx at z — 0, 
with its width set by the the observed scatter in this 
relation. 

For the highest values of psfr that are considered Fig- 
ure [2] (which would require significant incompleteness in 
current SFR measurements), not much evolution with 
redshift would be required to saturate the bound from the 
SXB. For the lower values of psFRj however, an order-of- 
magnitude increase in cx or more would be needed. At 
least for the galaxies with the highest SFRs, such signifi- 
cant evolution is marginally ruled out (jCowie etaDl20Tl 
Section |2~T|. 

3.2. Quasar Absorption Lines 

H I Lya forest and metal line absorption measure- 
ments from quasar sp ectra have been meas ure d to red- 
shifts as high as z w 7 (|Mortlock et al J 1201 It see lMeiksm] 
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Fig. 3. — Effect of hard sources on the EUV and X-ray back- 
grounds at z = 2.5, as well as constraints from the H I and 
He II Lyman forests a nd from the metal complexes studied in 
lAgafonova et all (2007). The fiducial background model discussed 
in the text in which quasars dominate the emissivity is given by 
the black solid curve in both panels. The other models in the top 
panel assume that an extra component of the emission was released 
at z = 10 and such that it saturates the 2 = unresolved SXB in- 
tensity with a = (dashed) and a = — 1 (dotted). Also shown 
is the upper bound for the z = unresolved SXB under the most 
conservative assumption that all contributing emissions originated 
from z > 2.5. The other models in the bottom panel assume that 
10% or 1% of hydrogen ionizing photons originate from a harder 
source than quasars. As in the top panel, the dashed curves as- 
sume a = and the dotted curves assume a = — 1 (shown only for 
the 10% case). Lastly, the ionization potentials of O VI and O VII 
are marked. O VI is the metal line that is sensitive to the highest 
energies for analyses at z = 2.5. 



120091 for a recent review of these diagnostics). These 
measurements are sensitive to the coeval soft X-ray back- 
ground. Unlike the z = unresolved SXB, constraints 
derived from high-redshift absorption lines avoid con- 
tamination from lower redshift, faint emission. How- 
ever, existing measurements of quasar absorption lines 
at z > 3 do not unambiguously constrain these radia- 
tion backgrounds as typically few lines are detected for 
each absorption-line complex (in part because the Lya 
forest is so saturated). Therefore, we concentrate on the 
constraints from z « 2.5 - the "sweet spot" for quasar 
absorption studies. 

H I and He II Lya forest measurements at z w 2.5 
constrain the the H I and He n photoionization rates, 
Thi and ThcII- We define r\ as the local ratio of the 
H I to He II column density, which is proportional to 
rm/rHcii in photoionization equilibrium. These quanti- 
ties can be related to the angular-averaged intensity per 
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where the second and fourth lines use c.g.s units, abk 
is the spectral index of the EUV background, and 
ax is the photoionization cross section at the ioniza- 
tion potential of ionic species X. Current measure- 
ments find r m ~ 0-5 - 1 x IP" 12 s" 1 and 77 ~ 80 
at z = 2.5 (IBolton et all 20051: iFaucher-Giguere et al.l 
2008bl iWorseck et al.l 120111 ). The error bars at 1 and 



4 Ry in Figure |3] show the constraints on J„ using equa- 
tions © and © with T H i = 0.75 (±0.25) x 10~ 12 s _1 . 

Metal absorption line measurements in quasar spec- 
tra also constrain the z ~ 2.5 EUV background (e.g., 
lAgafonova etafl [20071: iFechnerl f20lTT) . The red shaded 
trapezoids in Figure [3] are the best-fit EUV background 
spectrum for two z ~ 2.5 me t al ab sorption complexes 
presented in lAgafonova et al.l (|2007l) . where « 12 ab- 
sorption lines were simultaneously fit for each complex. 
We have specified the normalization of the shaded re- 
gion, which is unconstrained owing to degeneracy with 
the level of enrichment. In addition, the finite width of 
the shaded regions is for viewing ease and reflects a fac- 
tor of w 2 shift in this constraint. It does not reflect the 
uncertainty in this measurement a s error bars were not 
provided in lAgafonova et all (|2007t) . 

To investigate the impact of additional hard sources of 
emission on the EUV background and how much of such 
radiation the quoted absorption line constraints allow, 
we construct models for the 10 — 10 4 eV radiation back- 
ground. In these models, the mean specific intensity of 
the meta-galactic radia tion background at frequency 
and r edshift zq is (e.g.. IPeebleslll993T : IHaardt fc Madaul 
199G) 



(1 + z y 

47T 



dr]ae(is 1 z)e- T ^ ( » Z0 ' z \ 



(6) 

Here, dr\ = cdz/H(z) is the differential of the confor- 
mal distance, e{y, z) is the specific comoving emissivity 
at a frequency of v = vq(1 + z)/(l + zq), exp[— r e s] is 
the probability that a photon is transmitted as it travels 
between z and Zq, and 



Toff 



dz 



d/Vi 



HI 



dN 
dzdNm 
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(7) 



dAf / dzdNni is the H I column density distribution (con- 
strained by quasar absorption- line measurements), and 
t(Nbi) is the H I and He II continuum optical depths of 
a system with column ATjji at frequency v$ (l+z)/(l+z ). 
Equations © and (J7J) are exact in the limit of uncorre- 
lated (Poissonian) absorbers. 

We parameterize dN /dzdNm as a power-law w i th in- 
dex —1.5 (which is simplistic; iProchaska et al.l 120091 
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but our qualitative results are weakly sensitive to 
this choice). 9 In addition, we assume quasars with 
a spectral index of —1.7 dominate the emissivity, 
contributing a comoving emissivity at 1 Ry of 5 x 
10 24 erg s _1 Hz" 1 cMpc~ 3 at z = 3, and where this emis- 
sivity scales as (1 + z)~ 2 with redshift, in approximate 
agreement with estimates for th e contribution of AGN 
at z > 2.5 (|Hopkins et al.ll2007t ). (We ignore starlight, 
which may contribute as well; e.g. iFaucher-Giguere et al.l 
20083)- 

The top panel in Figure [3] shows the impact of hard 
sources on the EUV and X-ray backgrounds at z = 2.5. 
The black solid curve in both panels is the fiducial back- 
ground model that assumes only quasars. The other 
two background model curves in the top panel assume 
that additional emission was released at z = 10 with 
a = — 1 (dashed curve) and a = —1.5 (dotted curve), 
with its amplitude set such that it satur ates the unre- 
solved SXB. The red upper bound is the iMoretti et al.l 
(2003) constraint on the SXB (i.e., the z = unresolved 
0.5 - 2 keV SXB multiplied by [1 + z} 3 and blueshifted 
appropriately). The minimum energy that is impacted 
by intergalactic absorption, ~ 1 keV, is weakly model 
dependent, as it scales with the amount of hydrogen and 
helium in our models to the one-third power. Figure |3] 
shows that high-redshift emission at energies less than 
< 1 keV - energies at which absorption line measure- 
ments are most sensitive - is absorbed out by intervening 
gas by z = 2.5. Thus, the top panel in Figure [3] demon- 
strates that the z = unresolved SXB provides tighter 
constraints on relic radiation backgrounds from z ^ 2.5 
than z = 2.5 quasar absorption line measurements since 
it probes energies that are unabsorbed. 

The bottom panel in Figure [3] shows the J„ in the 
fiducial background model (solid curve) and for cases in 
which an additional 10%/1% of hydrogen ionizing pho- 
tons at 1 Ry originate from a hard source with a = 0/ — 1 
(dashed curves/dotted curve, shown only for a = 0). 
These cases are all clearly ruled out by the SXB mea- 
surement: The SXB plus absorption line measurements 
strongly limit any coeval hard emission. We have as- 
sumed in the bottom panel that the total emissivity of 
the hard sources falls off in the same manner as is as- 
sumed for the quasars, (1 + z)~ 2 . However, we find 
that the qualitative results are unchanged if their emis- 
sivity were instead constant with redshift such that the 
10%/ 1% contributions to J v at 1 Ry would only have 

9 We assume dM /dzdNm is normalized so that the the mean 
free path of 1 Ry photons is 300 cMpc at z = 3 and so that 
it decreases as (1 + z)~ 2 with increasing r edshif t, in agreement 
with measurements (Prochaska et al. 2009, 2010). Furthermore, 
we assume that the average diffuse He II fraction is given by 
tanh([2 — 3]/1.5), although the results depend weakly on the 
assumed He II ionization history because the frequency that has 
optical depth equal to 1 scales as the He II fraction to the one-third 
power. After He II reionization and in He II bubbles, we assume 
a 4 Ry photon has mean free path of 50 cMpc at z = 3, which is 
consistent with other e stimates (e.g., Fauchcr-Gigucrc et aLl l2009t 
IHaardt fc Madau|[20l2T) . 

Our simple background model does not include the emission 
from Type 2 quasars, which contribute 50-75% to the resolve d 
X-ray background at > 1 keV (see IHaardt fc Madau 2012). 
Accounting for unresolved Type 2 quasars would make it more 
difficult for these absorption line measurements to accommodate 
additional hard X-ray emission. 




redshift 



Fig. 4. — Exploration of whether quasars can ionize both the hy- 
drogen and helium and still be consistent with constraints on the 
cosmological reionization history. Shown are the average cosmo- 
logical ionization fraction in H II and He III for reionization with 
a quasar-like spectrum (a = —1.7), under different assumptions 
for the population's emissivity evolution and for the dumpiness of 
intergalactic gas. Here, eo = 8 X 10 — 25 erg s — 1 Hz -1 cMpc - 3 
(as motivated by Lya forest inferences) and Z3 = (1 + z)/4. The 
curves whose labels begin as "He III" are calculated with the same 
parameters as the "H 11" curves labelled directly above them but 
instead show the reionization history of He III. 

been satisfied at z = 2.5. 

In conclusion, the SXB places stronger constraints on 
hard high-redshift emission than z ^ 3 metal absorption 
lines under the assumption of power-law sources. Fore- 
ground absorption filters out high-redshift emissions at 
frequencies at which absorption measurements are the 
most sensitive. In conjunction, absorption diagnostics 
and the unresolved SXB show that a source population 
which contributes significantly to the ionization state of 
the intergalactic hydrogen and helium at z ~ 2.5 must 
also have a spectrum which cannot be much harder than 
that of low-redshift quasars. 



3.3. The Reionization History 

Constraints on hard high-redshift emission are also set 
by our knowledge of the cosmological reionization history 
of hydrogen and helium. The Lyman-a forest indicates 
that the hydrogen in the IGM was significan tly ionized 
by z ~ 6, when the Universe was 1 Gyr old ( [Fan et aLI 
|2006 | ; although, som ewhat lower redshifts are still possi- 
ble; [Mesingcr 2010). In addition, there is now significant 
evidence that He 11 reionization was ending at z = 2.7 
(Shull lFaT1l2Tjl0t iMcQuinnl [20091: iFurlanetto fc Dixon! 
2010: IWorseck et al. 1 1201 ID - well after hydro gen reion- 
ization and at the time the Universe was 2.5 Gyr old. 
We take the latter as a given for this discussion, noting 
that Gunn-Peterson troughs are se en in the He II Lya 
forest to redshifts as low as z = 2.7 (jShull et al.ll2010l ). It 
is very likely that thes e indicate large-scale He 11 regions 
(jMcQuinn et al.ll2009D . 

The Case A He ill — > He 11 recombination time is 
« 5.5 times shorter than for hydrogen at relevant tern- 



peratures and, at the cosmic mean density, is equal to 
0.7 times the age of the Universe at z ~ 3. Thus, if the 
sources that doubly ionized the helium shut off at a high 
rcdshift, the He in would have recombined into He II 
by z = 3 (except in deep voids; IVenkatesan et alj|2001t 
iMcQuinn fc Switzerll20l6T ). This shorter recombination 
time prohibits placing constraints from the observed end 
redshift of He II reionization on X-ray source models 
in which the sources shut off at sufficiently early times. 
It does, however, constrain models in which the sources 
were still active to lower redshifts. However, we argued in 
§3.21 that absorption line diagnostics plus the unresolved 
SXB exclude sources that are still active at z ~ 3 from 
having a harder spectrum than quasars. Thus, quasars 
are the most plausible sources this history can be used 
to constrain. 

We aim to calculate ionization histories for different 
source models. We define the variable C a to be the en- 
hancement in the recombination rate into ion a over a 
homogeneous Universe at a temperature of 10 4 K. At rel- 
evant redshifts, C a should increase with time (decreasing 
redshift) as structure forms and as the intensity of the 
ionizing background within ionized regions increases. 

Quasar absorption line measurements themselves 
place the empirical constraint Chi < 6 at z = 4 
(jMiralda-Escudl [20031: IMcQuinn et~aTI l201ll ) . In addi- 
tion, recent numerical studies ha ve found Chi < 5 in 
cosmological simulations at z < 6 ([McQuinn et al.|[2011[) 
and smaller val ues at higher redsh ifts, at least well af- 
ter reionization (jPawlik et all 120091 ). These studies also 
show that this quantity depends w eakly on the ioniz- 
ing background and also on redshift (|Pawlik et alll20Q{i 
IMcQuinn et al.ll2~011[) . To estimate Cneii numerically, we 
have p erformed a similar calculation to M cQuinn et all 
(|2011l ) in which EUV radiation is shone on all overdense 
regions within a cosmological simulation volume. These 
values yield ChoII = 1-5 — 3 over a large range of pho- 
toionizing backgrounds. We note that since dense regions 
tend to self-shield more easily at the ionization potential 
of He II than at this for H I, Cneii should always be less 
than Chi- 

Next, we take simple AGN emissivity histories that are 
either constant or power-laws in 1 + z. The emissivity 
of w I Ry photons has been measured at z ks 2 — 4 
and found to be just a few ionizing photons per hydro- 
gen atom per Gyr, not significantly in excess of what 
would have been required to reionize the hydrogen by 
z = 6 (iMiralda Escudd 120031 iBolt on fc Haehneltl l2007ri 
iFaucher-Giguere et al.ll2008al ). Thus, the comoving ion- 
izing emissivity cannot increase significantly from z ~ 10 
to z ~ 3 despite the massive growth of cosmological 
structures over this interval. 

Figure [4] shows the ionization histories for sources with 
a quasar-like spectral index of a = —1.7. These curves 
are calculated by solving for the ionization state of the 
hydrogen and helium given the attenuated incident spec- 
trum (tracked using many spectral bins). 10 The solid 
lines are the H II fraction as a function of redshift, and 
the other curves are the He in fraction. The different 

10 We set the electron density to be (1 — Yn c )p c (z)Q b /m p plus 
the electrons from ionic helium, which mimics an inhomogeneous 
reionization of hydrogen, where Yne ~ 0.24 is the He II fraction 
and pc(z) is the cosmic critical density. 



sets of curves take different clumping factors and emis- 
sivity histories, assuming that the sources turning on at 
z = 15 and where the emissivity is labeled in terms of 
e = 8 x 10~ 25 erg s _1 Hz " 1 cMpc~ 3 (roughly the value 
that is measured at z — 3; Haardt & Madau 2012). The 
red thick curves assume a constant emissivity of eo with 
redshift and Chi = Cneii = 2, consistent with recent nu- 
merical studies. In this case, the middle of reionization 
occurs at z = 9 (a bit on the low side of what i s pre- 
ferred by WMAP of 10.5 ± 1.2; ILarson et alllMlll) . H I 
reionization completes at z = 6, and He n reionization 
ends at z « 3.5. We require doubling the He n clump- 
ing factor - which is inconsistent with our estimates - 
for He II reionization to instead end at z = 2.7 (dot- 
ted curve). In addition, if the emissivity declines with 
redshift as (I + z)^ 1 from 2 eo a,t z = 3, this yields the 
blue curves, for which H I reionization ends at z = 6 and 
He II reionization ends at z = 4. These curves illustrate 
that the more likely case of a declining quasar emissiv- 
ity is more difficult to reconcile with He n reionization 
ending at z ~ 2.7. 

In conclusion, it is difficult for faint AGN (or, more 
generally, sources with a > —1.7) to be fully responsi- 
ble for reionization of both hydrogen and helium and to 
be consistent with the He II reionization history. Both 
hydrogen and helium reionization by faint AGN would 
require most of the following to hold: (1) that the AGN 
emissivity does not decline significantly with redshift to 
z ~ 10 in order to maximize the number of recombi- 
nations and, thus, the delay before He n reionization 
completes, (2) that the typical AGN spectrum does not 
harden with increasing redshift, (3) that the enhance- 
ment in the recombination rate is a factor of a few larger 
than in recent numerical studies, and (4) that hydrogen 
reionization ended at the lowest redshifts allowed by ob- 
servations, z«6. 

4. THE TEMPERATURE OF THE INTERGALACTIC 
MEDIUM 

This final section estimates the values for the in- 
tergalactic medium temperature that would have been 
achieved in plausible reionization models. Our goal is 
to understand whether an appreciable fraction of ion- 
izations by soft X-rays would have imprinted a distinct 
thermal signature. 

Photoionization heating is thought to be the dominant 
heating mechanism for the IGM at high redshifts (shock 
heati ng becomes important at z < 2; iCen fc Ostrikerl 
1999), and this heating depends on the sources' proper- 
ties, both their spectra and luminosities. Once a region 
was reionized, it would have been heated to a tempera- 
ture that depends on the incident ionizing spectrum. Its 
subsequent thermal history is determined by processes 
that are well understood: Compton cooling off the cos- 
mic microwave background, adiabatic heating and cool- 
ing, atomic cooling, and the photoionization heating of 
the residual neutral gas. 

The timescale for < 3 x 10 4 K ionized gas at 
the cosmic mean density to cool by a factor of 2 
at z < 10 is of the order of the Hubble time. 
Thus, many studies have noted the possibility of do- 
ing archeo logy via z = 2 — 6 temperature mea- 
surements (iMiralda-Escude fc R ecs 1994; iHui fc Gnedinl 
[19971 IHui fc Haimanl 120031: iTittlev fc Meiksinl 120071; 
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iFurlanetto fe Ohl H009). Roughly, temperatures of > 
10, 000 K at z — 6 indicate that part or all of the cos- 
mic gas was ionized within a redshift interval of Az » 3 
from this redshift because Compton cooling efficiently 
cools the gas that was ionized at z > 9 below such 
temperatures (|Hui fc Haimanl 120 03). Temperatures of 
T > 10,000 K at z « 3 likely indicate that helium 
rcionization occurred near this redshift as most of the 
thermal energy from hydrogen reioni zation has been lost 
(jTheuns et al.l [200l iHui fc Haimanl [200l . These con- 
clusions depend little on whether hydrogen reionization 
heated the gas to 15,000 K or 50,000 K. More pre- 
cise conclusions can be garnered with knowledge of the 
amount a region was heated to at reionization. 

There have been a few studies that investigated 
the impact of reioniza ti on on the IGM temperature. 
IMiralda-Escude fe Reed (|1994f) argued that if reioniza- 
tion occurred over a significant fraction of a Hubble 
time, collisional cooling would always result in signifi- 
cant cooling behind the hydrogen ionization front, lim- 
iting the temperature that the gas can attain. We at- 
tempt to make th is argu ment more quantitative here. 
ITittlev fc Meiksinl (|2007[ ) performed calculations that 
took into account the inhomogeneous nature of the IGM 
on the intergalactic temperature d uring reionization. 
Howev er, the planar source model in ITittlev fc Mciksin 
(|2007f ) enhances the prominence of radiative transfer ef- 
fects owing to inhomogeneities in the IGM and does not 
treat different front speeds (what we focus on here). 
IVenkatesan fc Benson! (|2011l ) investigated the tempera- 
ture profile in the H II region around a single galaxy, 
including stellar and black hole emissions, but investi- 
gated the heating just immediately around the galaxy 
and the did not investigate how much a typical parcel of 
gas was heated during reionization. Semi-analytic stud- 
ies of reionization have circumvented knowledge of how 
much heating occurs at the ionization front by leaving 
the post-front temper ature as a fixed parameter (e.g., 
IFurlanetto fc Ohll2009l there 30,000 K), and these stud- 
ies followed the subsequent evolution once a region was 
reionized. A few radiative transfer simulations of reion- 
ization on cosmological scales have a lso attempted to 
follow the heating from this process (jTrac et al.l 120081 : 
iFinlator et al.ll20lil) . but each study investigated a sin- 
gle spectral model and, since they did not resolve the 
ionization front, did not capture the amount of cooling 
at the front. No study has provided a systematic un- 
derstanding of how the thermal history depends on the 
spectrum of the sources and the morphology of reioniza- 
tion. 

We have developed a ID radiative transfer code that 
performs radiative transfer on the light-cone to calculate 
the ionization and temperature of the gas around a point 
source. This code uses the photoionization cross sections, 
recombination and col lisional ionizati o n co efficients, and 
cooling rates given in IHui fc Gnedhl (IT9971) . 11 We also 
include secondary heat ing and ionization processes, using 
the fitting tables from IFurlanetto fc Stoeverl (|2010D , and 
recombination radiation. In particular, radiation from 
recombinations are included in the approximate manner 

We do not employ crude power-law approximations for the 
cross sections, as is common in similar radiative transfer codes to 
speed up the computation. 
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Fig. 5. — Top panel: Temperature of the IGM around point 
sources with luminosities of 10 54 ionizing photons per second and 
with spectral indexes of a = —0.5, — 1.5 and —3. Each case is 
shown 10 Myr and 100 Myr after its source had turned on at z = 8 
in a fully neutral, homogeneous IGM. Bottom Panel: The same but 
where the emission is set to zero above 4 Ry so that He II cannot 
be reionized. In both panels, the vertical bars mark the location 
of the ionization front and are labeled by the velocity of the front. 
As the front's velocity slows, the amount of cooling increases and 
the difference in the temperature between the different spectral 
models decreases. The dotted curve in the top panel is the same 
as the a = —3 case after 100 Myr except there is a second, harder 
component with a = —0.5 and comprising 10% of the ionizing 
photons. A luminosity of 10 54 ionizing photons s —1 is a few order of 
magnitudes less than that of the brightest z ~ 6 quasars and could 
represent the emission from a faint AGN or from an accumulation 
of galaxies (see text). 



by adding these photons to the outward propagating ray. 
We have tested our code against the example ca lculations 
in the appendix of lBolton fc Haehneltl (|2007a|) . 12 

We ran ID radiative transfer calculations around a 
single point source in an initially neutral, homogeneous 
medium. We assume spherical symmetry, a luminosity 
of 10 ionizing photons per second, and that the source 
turned on at z = 8. We ignored adiabatic and Compton 
cooling as this prevents cooling well behind the front, 
aiding interpretation. We will argue that a luminosity of 
10 54 photons s _1 results in front speeds after 100 Myr 
with similar velocities to those expected during cosmo- 
logical reionization. However, our results can essentially 
be rescaled to other source luminosities by taking the 
temperature at the radius with the same incident flux, 
L/r 2 . 

Figure [5] shows the temperature profile from this calcu- 

12 These calculations must resolve the width of the H I ionization 
front or (cthi njj) m 0.7 Z~ 3 proper kpc, an impossible requirement 
for a non-adaptive 3D radiati ve transfer simulation on cos molog- 
ical scales (~ 10 Mpc). (See Cantalupo & Porciani 2011 for an 
adaptive 3D algorithm that can capture this physics.) In all of the 
calculations presented here, the front width is resolved with f» 10 
elements, which we find is sufficient for convergence. 
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lation for sources with spectral indexes of a = —3, —1.5 
and —0.5. The bottom panel is the same as the top panel 
except that a spectral cutoff is assumed above 4 Ry, such 
that the helium cannot become doubly ionized. Such 
a cutoff is meant to represent the case in which radi- 
ation from the starlight reionized the Universe. Com- 
monly, this case is assumed to have a very soft spec- 
trum with —a = 3 — 4, but recent spectral models are 
consistent with harder spec tra (e.g., see discussion in 
iFaucher-Giguere et al.1 12009). and absorption within the 
host-galaxy should also harden the emission that reaches 
the IGM. ' 

In both panels in Figure [SJ the thin dashed curves show 
the temperature after 10 Myr, when the ionization front 
is moving at 78 cMpc Gyr -1 , and the thick solid curves 
are the temperature after 100 Myr, when the front is 
moving at 17 cMpc Gyr -1 . The vertical bars mark the 
positions of the front (the H II region radius). The tem- 
perature immediately behind the ionization front is typ- 
ically larger at earlier times, when the front was moving 
faster. This temperature difference owes to the impor- 
tance of collisional cooling at the ionization front - the 
longer the front overlaps with a single gas element the 
more that element can cool (as the cooling rate peaks 
strongly at a hydrogen fraction of 0.5). Large front veloc- 
ities enable reionization to heat the IGM to > 30,000 K 
near bright sources (as can be seen in Fig. [5j smaller 
radii have larger temperatures). 

What speeds are anticipated for cosmological ioniza- 
tion fronts? It is thought that the Universe was ion- 
ized by small galaxies with a number density of n g ~ 
1 cMpc -3 , each producing on the order of 10 51 ion - 
izing photons per second (e.g., iMcQuinn et al.l 120071) . 
The front speed for an accumulation of galaxies with 
n g = 1 cMpc -3 in a ~ 5 cMpc H n region will be similar 
to that of the 100 Myr curves in Figure [SJ An alter- 
native estimate for the ionization front speeds is given 
by cosmological simulations of reionization. These sim- 
ulations show that reionization occurs over / ~ 0.5 of 
a Hubble time and with characteristic bubble sizes of 
Rb ~ 10 cMpc, implying front speeds of f~ 1 H~ 1 Rb ~ 
20 cMpc Gyr -1 . In detail, these simulations find that 
the fronts move slowly at first around individual sources 
but accelerate as more sources form and larger bubbles 
grow , having front speeds closer to 80 cMpc Gyr -1 by the 
end Jlliev et all [2001 IMcQuinn etUl [20071 iTrac fe Cenl 
|2007|) . Such acceleration would result in the regions that 
were ionized last being the warmest, especially for the 
harder spectra that we consider. However, the ioniza- 
tion front speeds at the end of reionization in these sim- 
ulations are almost certainly too fast, as these simula- 
tions are not capturing Lyman-limit systems, which ab- 
sorb ionizing photons and slow the growth of H II b ub- 
bles (iFurlanetto fc Ohl [20051 : IChoudhurv et al.l 120090 . In 
addition, the emissivity prescribed in these simulations 
tends t o overshoot the observed i onizing emissivity at 
z = 6 (jBolton fc Haehneltl l2007bf) , which also implies 
that reionization is happening too quickly in the sim- 
ulations. The longer the duration for reionization, the 
smaller the range of post-front temperatures. 

Thus, if the front speeds are of the order of ~ 
20 cMpc Gyr -1 , as we have argued, then the heating be- 
tween the different models considered in Figure [5] is not 



substantially different, ranging between 17 — 25 x 10 3 K 
for the wide range of spectra that we considered. Fac- 
tors of 4 faster speeds result in 2 — 3 x 10 4 K (thin 
dashed curves in the top panel of Fig. [5]), whereas slower 
speeds result in a slightly narrower range. Note that 
the temperatures for the a = —1.5 case can be larger 
than the harder a = —0.5 case. This difference owes 
to the larger amount of energy that was lost in sec- 
ondary ionizations/excitations in the a — —0.5 case, 
and it also owes to when the helium was reionized. If 
helium is doubly ionized ahead of the hydrogen ioniza- 
tion front (as for a = —0.5), then the imparted en- 
ergy is more significantly radiated away at the time the 
hydrogen ionization front pas ses by collisional cooling 
(jMiralda-Escude fc Rees1[l99l . 

Note that the case with reionization by sources with 
a = —0.5 is ruled out by the SXB f £|3.1|) . However, 
the scenario where ~ 10% of the ionizations occurred via 
such a hard background is still allowed. The dotted curve 
in the top panel is this case after 100 Myr where 90% of 
the photons are emitted with a = — 3 and 10% with 
a = —0.5. The temperature of the IGM is only slightly 
increased in this case over reionization by a source with 
a = -3. 

For reionization by stellar radiation - the most prob- 
able scenario -, it is unlikely that much > 4 Ry radi- 
ation would have been produced (and, if it were, that 
it would escape into the IGM). However, as mentioned 
above, there is significant uncertainty in the spectrum of 
the radiation that ultimately escapes a galaxy at 1—4 Ry. 
The bottom panel shows the temperatures in the case 
where the emergent spectrum has a = —0.5,-1.5, and 
—3 at 1 — 4 Ry. This panel demonstrates that, despite the 
uncertainty in the emergent spectrum for stellar sources, 
there is little uncertainty in the temperature to which 
the IGM would have been heated in this scenario. 

Our temperature profile calculations have been done in 
the limit of a homogeneous universe. We argue in what 
follows that density inhomogeneities should not signifi- 
cantly alter our results. In the limit of no recombinations 
(relevant for most inhomogeneities), the time an ioniza- 
tion front requires to pass a given point does not depend 
on these inhomogeneities. While in this limit the front- 
crossing time is the same, the cooling rate per particle is 
proportional to density so that overdense regions would 
have cooled to somewhat lower temperatures as the front 
crossed. In high enough density systems, recombina- 
tions can become important (which happens in overden- 

sities of ~ 100 rjj/ T 3 , in "L yman-limit systems" ; ISchavd 
120011: IMcQuinn et al]|2011f) . The continuum absorption 
of these systems can affect the post- front temperature, as 
they will slow the front (allowing for more cooling at the 
front) and harden the incident spectrum (so that more 
heat is injected per ionization). 13 Current thinking is 
that the number of ionizing photons required to ionized 
the un iverse is not much g reater than one per hy drogen 
atom ([Pawlik et al.l [20091: IMcQuinn etaLl f20Ll . This 
thinking is consistent with the finding that the Universe 

13 In all modern models for reionization many sources resided in 
an H II region during the bulk of reionization. In this picture, most 
sources' radiation incident on a gas parcel at the ionization front 
will miss the Lyman-limit absorber that obscures a given sightline, 
lessening the impact of such radiative transfer effects. 
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was not producing many more than one ionizing photon 
per Hubble time per baryon at z = 6 (iMiralda-Escudel 
12003 IBolton fc Haehnelti [20075 IMcQuinn et all 1201 ID . 
Thus, a large fraction of the emitting ionizing photons 
were unlikely to be absorbed in dense systems and change 
the character of the radiation field that reionized the 
IGM. However, even in a universe where there were a 
substantial number of absorptions in Lyman-limit sys- 
tems that hardened the incident spectrum, our calcula- 
tions for different spectral hardnesses show that the range 
of possible post-front temperatures is still small. This 
is particularly true for the (most plausible) case where 
reionization owes to < 4 Ry photons. 

Recent analyses have constrained the temperature that 
the IGM was he ated to by hydroge n reionization. The 
measurements of IBecker et all (|2011l ) constrain the IGM 
temperature to ~ 5 - 9 x 10 3 K at z » 4.5, and they are 
consistent with reioniz ation heating the IG M to 20, 000 K 
at z > 8. In addition. IBolton et al.l ()2011D measured the 
temperature in the proximity regions of seven quasars 
at Z Ri6, finding T = 12 - 20 x 10 4 K at 95% C.L. 
IBolton et al.l (|2011[ ) concluded that this temperature in- 
dicates that He n was ionized by the quas ar and that the 
hydro gen was reionized at z > 8 (see also Raskutti et al. 
120121 who accounted for more complex reionization his- 
tories) . Since we predict a small range in post-front tem- 
peratures with T « 20, 000K, these constraints are also 
consistent with all the models considered here - at least 
for the front speeds we argued are most applicable and 
as long as much of reio nization occurs at z > 8, the cur- 
rently favored picture (jLarson et al.ll20TTI ). 

5. CONCLUSIONS 

Several studies have suggested that the production of 
soft X-ray photons was more efficient per unit star for- 
mation rate at z ~ 10 than it is at z ~ and that these 
photons cont ributed s i gnificantly to the reionization of 
the U niverse (lOhll2001l:lMadau et al.l l2004l; iMir abel et al.l 
1201 It Uohnson fc Khochfarll2011l ). In addition, reioniza- 
tion by quasars may not yet be ruled out. While it 
has often been reasoned that quasars could not have 
reionized the Universe, this rationale has traditionally 
relied on observations of the decreas ing number density 
of bright, L > quasars at z > 3 (M adau et alj 11999 : 
iFaucher-Giguere et al.l [2008ah . This paper presented a 
census of constraints on high-redshift soft X-ray photon 
production. 

We discussed the constraints from (1) the unresolved 
soft X-ray background (SXB), (2) intergalactic metal line 
observations, (3) the late reionization of helium, and (4) 
the temperature of the intergalactic medium. We found 
that the unresolved SXB places interesting bounds on the 
contribution to reionization of sources with a spectral in- 
dex in intensity per unit frequency, a, greater than — 1. 
This bound limits the total number of ionizing photons 
radiated per hydrogen atom to be < 1 from several poten- 
tial sources: inverse Compton scattering off supernova- 
accelerated electrons, HMXBs (under empirically moti- 
vated spectral models), and dark matter annihilations. 
It limits the contribution to be < 1 in the latter two 
cases. All of these sources have been proposed as impor- 
tant contributors to hydrogen reionization. However, we 
showed that the SXB is consistent with quasars reion- 
izing the Universe, correcting a common misperception. 



In addition, we showed that if high-redshift galactic X- 
ray emissions account for the unresolved SXB intensity, 
significant evolution in hard X-ray luminosity-SFR rela- 
tionship from that observed at low redshifts is required 
(unless there was more than an order of magnitude addi- 
tional star formation at z ~ 8 than has been observed). 

We showed that z m 2.5 intergalactic metal line obser- 
vations (which constrain the coeval EUV background) 
are consistent with quasar- like sources with a s» —1.7 
dominating the extragalactic background. Metal absorp- 
tion line observations do not allow a harder component 
than quasars to source more than half of the ~ 200 eV 
background at z = 2.5. In addition, if a source pop- 
ulation with a = — 1 contributed even just 10% of the 
1 Ry emissivity at z — 2.5, it would overshoot measured 
bounds on the unresolved SXB. 

We found that the contribution of faint AGN to high- 
redshift ionizing backgrounds is most constrained by our 
knowledge of the reionization history of hydrogen and he- 
lium (but recent luminosity function measurements also 
excl ude much of the para meter space for AGN reioniza- 
tion; iWillott et afll2010bH . We showed that quasars can- 
not finish ionizing the helium until z = 2.7, while also 
ionizing the H I at z > 6, unless theoretical estimates 
for the dumpiness of intergalactic gas are low and the 
EUV emissivity of quasars/ AGN stays constant or in- 
creases with increasing redshift. 

Lastly, we investigated models for the heating of the 
IGM owing to reionization with a ID radiative transfer 
code. Our calculations account for the finite speed 
of light as well as secondary heating and excitation 
processes, effects that had previously either been ig- 
nored or treated in an approximate manner, but are 
in fact important to include when soft X-ray photons 
contribute to reionization. We argued that the bulk 
of the intergalactic gas would have been heated at 
reionization to temperatures greater than 17, 000 K 
and less than 25, 000 K for reionization by a wide 
range of spectral models (with an even narrower range 
if the sources cannot doubly ionize the helium). The 
narrowness of this range owes to the efficiency of cooling 
at the ionization front for plausible front speeds. The 
higher temperatures in this range did not necessarily 
imply that the gas was ionized by the hardest spectra, 
as our intermediate a = —1.5 case resulted in the 
highest post-front temperatures. In addition, the speed 
of the ionization front was as important for determining 
the post-front temperature as the emitted spectrum. 
Therefore, it would be difficult to distinguish models 
for the sources of reionization based on measurements 
of the intergalactic thermal history. Instead, the timing 
and duration of reionization is most important for estab- 
lishing the average temperature of the post-reionization 
IGM. 
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